clear 
use "/Users/marclerchenmueller/Dropbox/BMJ_Gender_Promotion/BMJ_Dataverse/HDV_data/clinical_articles_90pct_genderdesig.dta" 

***************POSITIVE FRAMING ANALYSES*******************************************************************************************************************************************************************

***Linear probability models for clinical research, overall and stratified by journal impact factor category 

areg i_positive i.l_female pct_female n_authors med_mesh_pctile/*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if i_clinical & pubtype == "Article", a(issn)
margins i.l_female
margins r.l_female

areg i_positive i.l_female pct_female n_authors med_mesh_pctile/*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if i_clinical & (jif > 10) & pubtype == "Article", a(issn)
margins i.l_female
margins r.l_female
**FIGURE 4

areg i_positive i.l_female pct_female n_authors med_mesh_pctile/*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if i_clinical & (jif <= 10) & pubtype == "Article", a(issn)
margins i.l_female
margins r.l_female
**FIGURE 4

*Including direct interaction

reg i_positive i.l_female##i.high_jif pct_female n_authors med_mesh_pctile/*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if i_clinical & pubtype == "Article", cluster(issn)
margins i.l_female if high_jif == 0
margins i.l_female if high_jif == 1

***************CITATION ANALYSES************************************************************************************************************************************************************************

**Note: When the dependent variable is log transformed, the coefficient on the dummy variable for positive framing is to be exponentiated to get to %change in citations as presented in article

areg log_cites i_positive i.l_female pct_female n_authors med_mesh_pctile /*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if (jif > 10) & (i_clinical) & pubtype == "Article", absorb(issn) 
estimates store clinical_high_jif
estimates table clinical_high_jif, eform  
**FIGURE 5

areg log_cites i_positive i.l_female pct_female n_authors med_mesh_pctile /*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if (jif <= 10) & (i_clinical) & pubtype == "Article", absorb(issn) 
estimates store clinical_low_jif
estimates table clinical_low_jif, eform  
**FIGURE 5

*Sensitivity analysis: Log transformed linear model with Scopus cites

areg log_scopus_cites i_positive i.l_female pct_female n_authors med_mesh_pctile /*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if (jif > 10) & (i_clinical) & pubtype == "Article", absorb(issn)
estimates store clinical_high_jif
estimates table clinical_high_jif, eform  

areg log_scopus_cites i_positive i.l_female pct_female n_authors med_mesh_pctile /*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if (jif <= 10) & (i_clinical) & pubtype == "Article", absorb(issn)
estimates store clinical_low_jif
estimates table clinical_low_jif, eform  

*Sensitivity analysis: Poisson regressions 

poisson cites i_positive i.l_female pct_female n_authors med_mesh_pctile /*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if (i_clinical) & pubtype == "Article" & jif > 10, irr vce(robust)

poisson cites i_positive i.l_female pct_female n_authors med_mesh_pctile /*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if (i_clinical) & pubtype == "Article" & jif <= 10, irr vce(robust)

poisson scopus_cites i_positive i.l_female pct_female n_authors med_mesh_pctile /*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if (i_clinical) & pubtype == "Article" & jif > 10, irr vce(robust)

poisson scopus_cites i_positive i.l_female pct_female n_authors med_mesh_pctile /*
*/ fe_2003 fe_2004 fe_2005 fe_2006 fe_2007 fe_2008 fe_2009 fe_2010 fe_2011 fe_2012 fe_2013 fe_2014 fe_2015 fe_2016 fe_2017 if (i_clinical) & pubtype == "Article" & jif <= 10, irr vce(robust)


